Correlations between Maxwell's multipoles for 
gaussian random functions on the sphere 
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■ Abstract. Maxwell's multipoles are a natural geometric characterisation of real 

, functions on the sphere (with fixed £). The correlations between multipoles for 

gaussian random functions are calculated, by mapping the spherical functions to 
random polynomials. In the limit of high £, the 2-point function tends to a form 
^Nl \ previously derived by Hannay in the analogous problem for the Majorana sphere. The 

t> ' application to the cosmic microwave background (CMB) is discussed. 

o 
o 

O PACS numbers: 02.30.Px, 05.45.Mt, 98.70.Vc 

A striking feature of randomness is the emergence of structure arising from 
independent processes. Here, I describe a simple, universal correlation structure 
associated with statistically isotropic random functions on the sphere: the correlations 
between Maxwell's multipoles. In the multipole representation, any real spherical 



function $(6*, </>), with fixed £, may be represented by i unit vectors Ui, . . . , U£, 
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>: *(M)= E °mCT0) = g(tti-v)-(ucV) 7 = = (i) 

•i-H . m =-e V x + y + z 

for a_ m = (— l) m a* m (ensuring reality of $), (x,y,z) restricted to the sphere and C 
a numerical constant. The i vectors Ui are called Maxwell's multipoles |3J, and are 
determined uniquely (up to sign) [21 El 13- The significance of the Maxwell representation 
is that the multipoles rotate directly with the function, and are defined without reference 
to an external reference frame (unlike spherical harmonics). An arbitrary spherical 
function with fixed £, with its corresponding multipole directions, is shown in figure d 
The quantity discussed in this note is the correlation of the Maxwell multipole 
directions for completely random functions of the form (JTJ. Complete randomness means 
that the probability density function of $ depends only on the sum Ylm=o \ a m\ 2 ', the 
statistical distribution of the multipoles (which is automatically rotationally symmetric) 
is determined absolutely by this restriction. Here a gaussian distribution of the 
sum is chosen for calculations, which implies that the a m coefficients can be taken 
as independent (for m > 0), identically distributed gaussians (by analogy with the 
derivation of the Maxwell distribution in kinetic theory). 

| Present address: School of Mathematics, University of Southampton, Highfield, Southampton S017 
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Figure 1. A sample spherical function $ for I = 10. The radius of the plot is 
determined by &{9, <fi) plus a constant, and shaded according to the value of The 
Maxwell multipole directions for the function are also represented. 



The statistical morphology of cosmic microwave background (CMB) radiation is 
still little understood, and it has been suggested that correlations between Maxwell 
multipoles might reveal hidden structure in the data jSJ Ej. Here, I will derive 
the universal correlation function of Maxwell multipoles for a statistically isotropic 
spherical function, using independent gaussian random a m coefficients. The present 
analytic approach contrasts with the numerical Monte Carlo computations for the same 
isotropic gaussian distribution studied in the papers referred to above. 

The main tool used in the calculation is the Majorana representation of spherical 
functions (also called the 'stellar representation') [H IH1 El CHI HI] • In this representation, 
the function 0) with the form is represented as a polynomial in the complex 
variable £, 

/(o=£(-ir Lfj a ^ e+m ( 2 ) 

m=-i \ J 

The relationship with the sphere is realised by taking ( = exp(i0) tan(6 l /2) as a 
stereographic coordinate in the complex plane. Since each a_ m = (— l) m a* m1 the roots 
of the polynomial / occur in antipodal pairs Q,—l/Q; these are the stereographic 
projections of the Maxwell multipole vectors (whose sign is undefined) jl]. For complex 
functions on the sphere (with independent a m ), spinor-like Majorana directions arise 
from the roots of (J2J). (In the language of group theory, basis functions of the 
unitary (rotation) group SU(2) are mapped to those of the linear group SL(2,C) in 
this representation.) 

The identification of the spherical function with gaussian random coefficients, 
with the polynomial 0, reduces the problem to a calculation of the statistics of 
roots of random polynomials. The roots of spin-like random polynomials (SU(2) 
polynomials) have been studied, particularly in connection with quantum chaotic 
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systems [111112111311111131131131111121, with the distribution ©, but different 
restrictions on the identically distributed gaussian coefficients a m . 

The present calculation is similar to that of Hannay [H], who derived the Ap- 
point correlation function p k between the Majorana spinor directions for complex 
functions analogous to $ in equation (JTJ), that is, with no additional requirements on 
the coefficients a m . The 2-point function p 2 , in the large I limit, was found to have a 
simple universal form. Following the same method, Maxwell's multipoles are found here 
to limit to the same function. Readers not interested in the details of the calculation 
may skip to equation (0) and the following discussion. 

As stated above, $ in equation ((H) is assumed to have coefficients a m independent 
and identically gaussian distributed (which implies, but is not implied by, isotropy): 

(a m a n ) = 5 mn , (a m a n ) = ( 1) 5 m _ n . (3) 

All statistical information about the multipoles is in the zeros of the Majorana 
polynomial (J2J), which are accessed by averaging /j = f(Q) at the point Q, and its 
derivative f[ = df(()/d(\^ = ^. The joint probability density function P for k points 
fx, • • • , fk, fi, ■ ■ ■ , fh, fx, ■ ■ ■ , f' k , fx, ■ ■ ■ , fk is g iven b y the gaussian distribution 

p/ r f f* f* f' f' f'* f'*\ — ex P( — F • M • F/2) 
r Ui, • • • 5 Jk, J 1 , • • • 5 Jk, Ji) • • • > Jk, Jx i ■ ■ ■ > Jk ) ~ 7r k y/ det M ' 

where F = (fx, ■ . . , f k , fx, ■ ■ ■ , fk, fx, ■ ■ ■ , f k , fx, ■■■,/£), and M is the 4k x 4k hermitian 
correlation matrix M^- = (F*Fj). Following [llj, submatrices (of dimension 2k) of M 
will be identified: 

-(etc)- <0 

Nonvanishing terms (and their conjugates) appearing in M are JT] 

(fifj) = (i + (*Q) 2i , (fifj) = (C* - Q 2e , 

(ftf) = 2^0(1 + QQ 2t ~\ (tifj) = 2i(d - Q 2t -\ 

(ft/-) = 2C(i + QQ) 2e ~\ (/*/;> = -nd - q 21 -\ 

(ft ft) = 21(1 + 2££( j )(l + QQ 2l -\ (ffi) = -2t(2t - 1)(Q - Q 21 ' 2 , (6) 

following from averaging the complex polynomials (J2J) using ©• 

The fc-point function Pk(Cx, ■ ■ ■ , Cfc) is found using standard methods for finding the 
correlations of zeros of gaussian random functions (e.g. 20 ): 

Pfc (Ci, . . . , CO = WiMffl ■ ■■s(f k )5(f* k )\f[ ■ ■ ■ f k \ 2 ) 

= —7^= / d 4k F5(f 1 ) ■ ■ .5(fZ)\f ■ ■ ■ f k \ 2 exp(-F* ■ M _1 ■ F/2) 
7r Vdet M J 

™/ W|/! '"' il ' B ' ( - r ' r '- f ' /2 »' (7) 

where F' = (f[, . . . , f k , f[*, . . . f k ), and N _1 is the submatrix of M _1 in the position 
of C in M. By Jacobi's determinant theorem HH HI], N = C - B+A X B, and 
det M = det A det N. 
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The final part of the calculation is a gaussian integration by parts (i.e. the finite 
dimensional analogue of Wick's theorem), and the final result is (7T k V det A) -1 times a 
combinatorial term. This term is the sum over products of elements of N, one for each 
pairing of the 2k terms /{, . . . The elements of N in the product are those which 
correspond to the appropriate components of F', F'* in the quadratic form. 

The simplest nontrivial case is the 2-point function p%, which is 

p 2 = {N U N 22 + N 12 N 21 + AT 14 Ar 41 )/(vrVdetA) (8) 

It should be noted that, since the /j are correlated with each other (rather than 
simply with their conjugates), this combinatorial term is rather more complicated 
than the simple permanent found in |TI]; furthermore, the choice of some elements 
of N is not unique (for example, in the third summand in the numerator in (JHJ), 
N 14 = N 23 ,N 41 = N 32 ). 

The 1-point function p\ is simply the density of roots in the complex plane (counting 
each multipole direction, undetermined in sign, twice). On the sphere, this is uniform 
with value 2£/4n. The 2-point function p 2 can be simplified by taking the two points 
to be and r (with r = tan(0/2) for points with angular separation 6); no generality is 
lost since the distribution is rotationally symmetric on the sphere. From (|SJ), this is 

p 2 (0, r) = (^D^y'ifieD - Abuv - (b 2 + v 2 )(a - 1 - u 2 )} 

x[dD - 2cuv(a + 1 - u 2 ) - (c 2 + av 2 )(a - 1 - u 2 )} 

+ [2£D - 2cuv - buv(a + 1 - u 2 ) - v 2 (a - 1 + u 2 ) - bc(a - 1 - u 2 )} 2 

+ [wD - 2bcu - uv 2 (a + 1 - u 2 ) - bv(a - 1 + u 2 ) - cv(a - 1 - u 2 )} 2 } (9) 

with D = det A = (a — 1 — u 2 — 2u)(a — 1 — u 2 + 2u) and, extending the notation of [11J, 

a = (1 + r 2 ) 2£ , b = 2£r, c = 2£r(l + r 2 ) 2l ~\ d = 21(1 + 2£r 2 )(l + r 2 ) 2l ~ 2 , 

u = r 2 \ v = -2£r 2£ -\ w = -2£{2£ - l)r 2 ^ 2 . (10) 

To find the 2-point multipole correlation function function on the sphere p S phere(#), for 
points with angular separation 9, this must be multiplied by the stereographic jacobian 
(1 + r 2 ) 2 /16 and substituting r = tan(#/2) [TT]. p sp h e re(#) is symmetric about ir/2, by 
antipodality of the multipole directions. The function may be normalised (i.e. is unity 
when points are uncorrelated) by dividing by the square of the density 2£j\. r n. When 
£ = 2, p S phere(#) can be found from equation (jHJ), 

agreeing up to a numerical constant with that found by [Zj. p S phere is plotted against 
9 in figure El for some choices of £. There is good agreement with the corresponding 
numerical results of 7\ (who plot against cos 9) . 

In [TT], it was found, for a general complex isotropic Majorana polynomial, in the 
limit of large £, that p 2 (0,r) tends to (£/2ii) 2 g(\/Jr) ) where 

g{R) = ((sinh 2 R 2 + R 4 ) cosh R 2 - 2R 2 sinh R 2 ) / sinh 3 R 2 . (12) 
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Figure 2. The 2-point correlation function on the sphere p S phere, plotted against 9 
(in degrees) for multipoles with separation 9 for (a) I = 3; (b) t = 5; (c) I = 10; (d) 
i = 100. The function has been divided by the square of multipole density £/2ir (so 
approaches 1 for large £), and the behaviour for small 9 is quadratic. 

In equation (jUJ), on writing r = R/ \/J and taking the large £ limit (in which u, v, and w 
approach 0), equation © factorizes to give the 2-point function of this general complex 
case (equation (11) in JI]). In the complex plane, this large £ limit of the general 
random Majorana polynomials is the universal 'chaotic analytic function' [To^ I19j. 

Different symmetries for the random polynomial coefficients have been investigated, 
such as for a m real [13]. In this case, the distribution of roots in the complex 
plane behaves singularly on the real axis (where roots appear in complex conjugate 
pairs), and approaches the chaotic analytic function behaviour far from these lines 
(where there is little influence from the conjugate, paired zero). (The generalization of 
random polynomials with coefficients of fixed argument ^H] have similar behaviour.) 
The distribution of roots of the Maxwell multipole polynomials is uniform on the 
(stereographic) sphere, and the difference with the chaotic analytic function polynomials 
is the presence of the antipodal root, whose influence diminishes as £ increases. 
Therefore, for high £, the correlations between Maxwell multipoles also limit to the 
form (fT2"|) . apparently by a power law. 

g(R) exhibits quadratic repulsion at the origin with a small maximum of 1.0531 at 
R = 1.4985, and thereafter approaches 1 exponentially. For high £, the maximum of p2 
(visible in figure EJ) is approximately at 9 = 3/2\/i. In this limit, the 2-point function 
of Maxwell multipoles is therefore constant (there is no correlation) apart from a tiny 
repulsion region near each multipole. This is analogous to the closely related system of 
a one-component plasma on the sphere [21], whose 2-point correlation function, in the 
thermodynamic limit, approaches a constant apart from repulsion very near the origin. 

Higher correlation functions between more Maxwell multipole directions, which are 
universal for completely random spherical functions, can also be calculated using the 
methods described here. For large £, these will tend to the same behaviour as high £ 
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random Majorana polynomials |!Hll9j. However, it is possible that lower i correlations 
might reveal more subtle structure, which may be compared with CMB data. 
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